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Abstract 



In this paper we investigate the possibility of finding models that, in the global 
minimum, break SUSY by a combined F and D-term effect. We find that if the su- 
perpotential is a cubic polynomial in the fields and the Kahler potential is canonical 
this does not happen. 

We then give two examples where the previous conditions do not apply and 
where the vacuum of the model breaks SUSY by combined F and D-term vevs: 
SQCD+singlet with Nq = Np + 1, and a baryon deformed ISS model. 

We briefly discuss the phenomenological viability of the models. 



1 Introduction 



There has been, after the construction of Intriligator, Seiberg and Shih [T] , an interest 
in studying models that break SUSY on a metastable ground state. Many of these con- 
sist in deformations of the original ISS model in such a way that R-symmetry is either 
approximate [2], [3] or spontaneously brokenjH [5J Ej. The breaking of R-symmetry is im- 
portant if upon integration of the messenger fields we want to generate Majorana masses 
for the supersymmetric standard model (SSM) gauginos. 

In the IR, most of these models can be described by some O'Rafeartaigh-like dynam- 
ics. Supersymmetry is broken because not all the F-term equations can be solved, often 
appealing to the rank condition. In many of these, unfortunately, gauginos were pre- 
dicted to be much lighter than scalars. This is not so much related to how one breaks 
R-symmetry, but because we are expanding the potential around its tree level minimum 

m- 

The reason why it's not good to have a large hierarchy between these masses has to 
do with the amount of fine-tuning one faces when trying to include a mechanism that 
generates \x and terms to achieve electroweak symmetry breaking. 

The understanding of why the ratio of gaugino and scalar masses is large in these 
models already gives some hints of ways to avoid it. One possibility [8j is to use the 
ISS construction, but expand the theory around a different vacuum which has a higher 
energy. One has to ensure that the vacuum is stabilized somehow, but the argument of 
the vanishing of gaugino masses does not apply. A similar possibility is of a hybrid model 
that mixes explicit and direct mediation in such a way that one can tune the parameters 
of the theory to decide how metastable the vacuum is [9] and change the gaugino/scalar 
mass ratio. This approach assumes a different kind of metastability than the one of the 
original ISS proposal, as in these models there is a SUSY vacuum that can be found 
within the limits of validity of the effective low energy theory. 

Another hypothesis, which we will follow here, is to consider the effects of switching 
on the gauge coupling and look for a combined effect of F and D-term SUSY breaking. 
Examples were found in the context of semi-direct mediation in [TO], [11] , or [T2l [13] . 

These theories of semi-direct mediation also have problems in generating large gaugino 
masses when compared with scalar masses. An approach to solve this problem has been 
given in|14j. 

One could also try to look for pure D-term SUSY breaking by introducing a Fayet- 
Iliopolous term. However [T5| this is not preferred. This is because such a term is not 
invariant under complexified gauge transformations. Under these, the F.I. term induces 
a Kahler transformation. This is a kind of transformation that has no impact on the 
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lagrangian of the system, but makes the supercurrent multiplet not gauge invariant. So, 
if such a term is not present in the high energy theory, gauge invariance of the supercurrent 
protects it from being generated at low energies. If it is there at high energies, this theory 
can only be coupled to supergravity if there is an exact global R-symmetry(in the SUGRA 
theory). It is, however, problematic to have a theory of gravity with a global symmetry, 
so F.I. terms should vanish. This argument does not apply to "field dependent F.I. terms" 
coming from "anomalous" U(l)'s p3], but only true (field independent) F.I. terms. 

The other question that arises is the why should the universe be in a metastable 
vacuum. One possible reason for this is that because of finite temperature corrections to 
the effective potential in the early universe would make this state lower in energy and 
thus the true vacuum of the theory [T71 [T8] . 

In this paper we will look at the possibility of having models of combined F and D 
term SUSY breaking on the global minimum of the (effective) theory. In section 2 we will 
give some constraints on classes of models that can break SUSY in this way. In Section 
3 we will study the case of SQCD with Nc > Np and a singlet, to see that the vacuum 
breaks SUSY with both F and D-terms. This example is, in some sense, complementary 
to [191 [20] that study SUSY breaking in models with vector-like matter and evade the 
Witten index argument. However, the minimum we find does not have any flat direction. 
In section 4 we will study a baryon deformed ISS model and see that the stabilization 
of the runaway by the Coleman- Weinberg contribution generates small D-terms; we also 
comment on the phenomenological viability of this model. In section 5 we present our 
conclusions. 



2 Some constraints on models that can have com- 
bined F and D-term SUSY breaking 

In this section we will show that if the superpotential is renormalizable and the Kahler 
potential is canonical, the model can only break SUSY through F-terms (globally), i.e. 
when minimizing the function Vd + Vf one always finds Vd = 0. To do this we will find 
a set of solutions that minimizes Vf and then show that it's always possible to chose one 
that solves Vd = 0. Since Vd > this is a global minimum of the model. 

Since we want to minimize Vf, we start by reviewing some known results about 
O'Raifeartaigh (O'R.) mo dels [ZD, [7]: 

Let us consider a renormalizable Wess-Zumino model with superpotential W (i.e. W 
is a degree 3 polynomial in the fields), the Kahler potential is canonical. Then, Vp = 
Y2i l^il 2 ; where Wi is the derivative of the superpotential with respect to the field 0* and 
indices are raised and lowered by complex conjugation. 
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Since we are assuming that the minimum of the potential breaks SUSY, not all the 
Wi can be made to vanish simultaneously. At the minimum, the gradient of the potential 
must vanish, WijW^ = (i.e. we are not considering cases where there are runaway 
directions, even though they might be interesting). 

The tree-level boson mass is: 

2 _ / M* F M F T* 
m z-\ T M F M* F 

where Tij = WijkW k , and M.p = Wij. In a consistent vacuum M\ must be positive 
semi-definite. 

We will now see that in the direction (W l , Wi) this scalar mass matrix has a zero 
eigenvalue: 

Since this mass matrix is positive semi-definite, this must vanish, otherwise it could be 
made negative by rotating the phase of W % . So W l J r i jW^ = 0. It's possible to show, by 
performing an expansion of the potential to 3rd order, and using this equation, that one 
actually has the stronger result: TijW^ = 0. 

If the superpotential is renormalizable, then the model has a classically flat direction 
parametrized by <p % = + zW i : 

5Wi = W ijZ W j + W ijk (zW j )(zW k ) = (3) 

This expansion is exact since higher derivatives of W vanish. So, the F-terms are constant 
along this direction. 

We will now consider the two possible situations: 

1. The F-terms do not break gauge symmetry; 

2. The F-terms break gauge symmetry; 

The first option is, for example, the case of the ISS model, the second case looks like 
the 3-2 or 4-1 models, except that in those cases the superpotentials are not 3rd order 
polynomials. 

In the first case the non- vanishing F-terms are all gauge invariant. This means that 
if we perform a complex gauge transformation on the vevs of the fields, the potential 
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remains invariant. In this case the symmetry group under which the potential Vp is 
invariant is enlarged from SU(N) to SL(N). This symmetry is enough to solve all the 
D-term equations. 

One should note that this is not a true complex gauge transformation as only the 
vevs of the fields, and not the fields themselves, are being rescaled. Meaning, if we were 
performing a complexified gauge transformation, invariance of the potential would be 
guaranteed by a shift in the lowest component of the vector superfield. In the case we 
are considering, the SL(N) symmetry of the vacuum would still be there if there was no 
vector superfield (i.e. gauge symmetry). 

We consider two cases to illustrate the point: one explicit example and a general 
argument. 

Take a ISS model where we have gauged a U(l) baryon symmetry. The model has 
gauge group U(l) and the flavour group is SU(6) (for this number of flavours the magnetic 
gauge group of the ISS is empty). The field content is: quarks, 0, that transform under 
the fundamental of color and flavour, antiquarks, 0, that transform as anti-fundamentals 
of color and flavour, and some mesons, $, that are color singlets and transform under the 
adjoint plus a singlet of flavour. 

The superpotential is: 

W = h(Tr[<j>4&] - n 2 Tr[<&}) (4) 
SUSY is broken by the rank condition: 

= Kh^i - ffiv) ( 5 ) 

Since the number of colours is less than the number of flavours, it's not possible to solve 
all these equations and SUSY is broken. The scalar potential is then minimized to be: 

V F = 5h 2 fi 4 (6) 

The moduli space up to global symmetries is given by: 

*=(o ° x ),0=(^),0=(^),with O O = / / (7) 

The D-term potential in a random point of the pseudomoduli space is 

Vd = y(l0o| 2 -|0o| 2 ) 2 (8) 
Under a complexified gauge transformation of the vevs, 
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And V D (a) = (e 2a |0 o | 2 - e^ 2Q |0 o | 2 ) 2 - This vanishes when a = ±Log[^-. Note that 
on the F-term moduli space, neither <fi or <fi can vanish. 

For the second example we just recall the argument given in[22] and note that it still 
applies if the F-terms are non- vanishing, but gauge invariant. 

For this particular example SU(N) is semi-simple, so it has a Cartan sub-algebra. The 
D-terms transform in the adjoint representation of the group, and it is always possible to 
find a group transformation that rotates any initial D-term into the direction given by 
one particular element ( call it T) of the Cartan subalgebra. This gauge transformation 
gives another (gauge equivalent) solution to the minimization of the potential Vp. 

This generator T = diag({ni}) (where diag({fii}) is a diagonal matrix with eigenvalues 
{/ij}), and 

D = ^.diag({fXi})4 (10) 

Where <p and dot product is used as a notation for the sum over all the fields with the 
generator in the appropriate representation. 

Because the F-terms are gauge invariant, we can actually complexify the parameter 
that defines the gauge transformations, and still have a (physically distinct) solution of 
the minimization of Vp. This argument only differs from [22] on the assumption that all 
the F-terms vanish. 

We now distinguish two cases: 

• all fii have the same sign. 

• the ^ take both signs 

Under a complexified gauge transformation the D-term changes as 

D = fafne 2 *^ (11) 

In the first case we just need to take 77 — > — 00 if all the \ii > or 77 — ► +00 if \ii < 0, 
to have a solution to Vd — 0. 

In the second case we note that the D-term can be written as a total derivative 

D = —(ke 2 ^ (12) 
2 or] 

The function tpie 2 ^ 1 '<fr l goes to +00 as 77 — > ±00 and so it has a minimum where the it's 
derivative must vanish. In this minimum the D-term vanishes. As known, this is why in 
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U(l) theories without an F.I. term or non-Abelian theories, pure D-term SUSY breaking 
is not possible. 

We now look at models where the F-terms are not gauge invariant. In this case the 
vacuum manifold does not have the enhanced symmetry, and the previous argument does 
not hold. However, the vacuum manifold contains at least the direction parameterized by 

ft = (0*)(°) + Z W\ 

Since not all F-terms are gauge invariant, the heavy vector superfield will not be SUSY; 
for example there will be mixings between gauginos and chiral fermion mass matrices. 

One can still investigate how different matrix elements are constrained by gauge in- 
var iance: 

Wi{D a Y = Wi(T a ))<p j = (13) 

Because the superpotential must be gauge invariant, note that this is an identity and is 
not valid only at the minimum of Vp. Because of this we can get further identities from 
differentiating the previous expression: 

w ij {D a y + w i {D a )i = o 
w lJk {D a y + w i3 {D a Y k + w lk {D a y; = Q 1 } 

At the minimum of the Vp, the previous expressions implies that: 

Wi(T a y k W k = -WikW^DY = (15) 

We will now see that the D-terms are constant along the pseudo-moduli direction 
fa = 0f) + z Wi. 

D a (z) = D a (0) + 2Re(z((f)f ) )(T a ) i j W j ) + \z\ 2 W i {T a )^W j = D a (0) (16) 
The linear term in z vanishes because off [TBI) , and the quadratic term in z vanishes by 

We will now show that even though the D-terms are constant along that particular 
pseudo-moduli direction, they can be made to vanish arbitrarily close to it. We will give 
first a specific example and discuss why this happens and then give a general argument. 

Let us take a model considered in [23]: we have two chiral/antichiral fields that trans- 
form under a U(l), and a singlet. The superpotential is: 

W = hS(4>X - yu 2 ) + mi 00 + m 2 xx (17) 
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Except for S, the untilded fields have charge +1 under the U(l) and the tilded fields 
have charge -1. The model has an R-symmetry with charges R(S) = 2, R{4>) = R(x) = 
0,R(4>) = R(x) = 2, and the superpotential is general. 

The F-term potential is given by 



V F = h 2 \(j)x - M 2 | 2 + m\\x\ 2 + m 2 2 \(f)\ 2 + \hSx + W| 2 + \hS<f> + m 2X \ 2 



The model breaks SUSY at tree-level, and the minimum is: 



V F = h 2 f i 4 if /i 2 <^ 

V F = 2m 1 m 2 /i 2 - ^ if /i 2 > n^m 



(19) 

— ii p ~2hT 



This as long as both m\ and m 2 are positive. This model has pure F-term SUSY breaking: 
In the first case no charged field gets a vev. The U(l) is unbroken and the D-terms vanish. 

In the second case the U(l) is broken. The vevs are: 

, y/m.2yj m 1 (h 2 fj, 2 — rn 1 rri2) 

t hmi 



7 \l m\(h 2 y? —ra\m,2)S 
<P = 7= 

-y/ mi(/i 2 /i 2 — mim2)S 
^* mi v 'm2 

_ \ /mi (h 2 fj, 2 — m\rri2) 
Y — 



(20) 



for any value of S. 

The flat direction is parameterized by: 



S = - Z 2apa 

/i 

y/m^^J m 1 {h 2 y 2 -mim,2) s 
ft ^ 
v, _ v (0) _ y^Vmti/iV-mimi) 
A A /, * 



If we replace these vevs into Vn, we get that Vb = hg 2 ^™ 1 ^L^ 1 "^ 2 h M ^ , which is 
generally non-zero unless mi = m-i- Note that this is a constant and independent of z, as 
we have seen it should be. 
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We will now deform the vevs by e in the following way: 

, _ yJ m 2\J mi(h 2 fi 2 -mim 2 ) 

t hmi 

7 \J mi(h 2 /x 2 -mim 2 )5 

mi Wrri2 



-^/mi(h 2 fi 2 -mim2)S ^ ^ (22) 



_ \/mi (h 2 fj, 2 — m\rri2) 
Y — 



In this vacuum 



„ (m 2 — m 2 .)^/ mi(fi 2 /i 2 -mim2) 



Vf = 2m 1 m 2 /i 2 - + (mj + m 2 )e 2 (23) 

AT 

V D = (24) 
So, if e — > 0, we solve Vd = and minimize Vp. 

More generically will show that one can find a small deformation e of the pseudo- 
moduli space, where the change in the D-terms is proportional to ez, while the change in 
the F-terms will only be proportional to e. One can then take e — > with ze fixed, to 
solve Vd = 0, and Vp will approach it's minimum value. This means that if we ignore 
D-terms the model breaks SUSY at tree-level through F-terms, but as soon as we include 
D-terms they either vanish identically on the pseudo-moduli space or there is a runaway 
to broken SUSY. 

The stabilization of this runaway direction has to be checked because in these models 
one has the extra contribution from gauge fields to the Coleman- Weinberg potential, and 
these are not positive definite. This fact can be seen, for example, in [23] and in the 
last example model we show. Some cases where some metastable vacua with both F and 
D-terms and no flat directions have been built 



We will consider that a gauge rotation of the fields has been made in such a way that 
there is only one non- vanishing D-term in a direction defined by an element of the Cartan 
subalgebra of the group, very much like when we dealed with the case that all the F-terms 
were gauge invariant. 

We will chose that the deformation is a complexified gauge transformation in the 
direction defined by the same element of the Cartan subalgebra. This transformation will 
generically change both V F and Vd- 

Let T be the element of the Cartan subalgebra, then T = diag({fii}). Under this 
complexified gauge transformation The F-terms change as 

Wi{z) = W j (z)(e-°* Tb )i = (W(0) t .(e- afl, ) i ) = (W t .dm^({e^})) i (25) 
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Where, again, we have used (p and dot product as a notation for the sum over all the 
fields and the generator is in the appropriate representation. 

Then, the change in the potential is 

V F {a) = Wldiag({e a>ii }).W « V F (0) + a?W* .diag({$}) W + 0(a 3 ) (26) 

This is because, as we've seen, the terms linear in the generator of the gauge group vanish. 
The previous example where we had a U(l) gauge group agrees with this expression. 

We now use the expression (fT2"|) evaluated at (</>*)' = <p % + zW l . 

We get that: 

D{z) = ^.diag({tH}).<f> (27) 

- 2/m(a)(0 t .diag({tf}).(j) + 2Re(z$ .diag ({//■}). W) + \z\ 2 W ] ' .diag({tf}) W) 

(28) 

When z is very large this is solved by: 

[a) 2\z\*Wi.diag({ tJ *}).W 1 j 

The solution a is small and the corrections in a 2 we ignored are negligible if we take 
z — > oo. Note that Wi^L 2 W % = ^ 4 |/ijH^*| 2 and should be greater than zero if SUSY is 
broken and the F-terms are not gauge invariant. 

The conclusion is that in this class of models it's not possible to have a model where 
the global minimum has combined F and D-term SUSY breaking. 

One should, of course, be careful to check if the description being used is correct and 
one doesn't go to a strong coupling regime or if one is using an effective theory, go outside 
its regime of validity. 



3 Adding non polynomial terms to the superpoten- 
tial 

In the last case, the discussion relied in the existence of the flat direction, so one question 
is: what happens if we do not have a flat direction. One such possibility is to consider 
the present model, but instead of a U(l) gauge symmetry, one considers a SU(N C ) gauge 
symmetry with iV c > 2. In this case we have a SQCD theory with 2 flavours and a singlet 
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and N c > Np. This theory has a dynamically generated term in the superpotential that 
comes from instanton contributions, and is not a polynomial of degree 3 (for simplicity 
we take N c = 3). 

The full superpotential is: 



W = hS Ox - /i 2 ) + mx00 + m 2 xx + —~ — : ~— (30) 

<?4x-x - x4<?-x 

At tree-level the superpotential is general and has an R-symmetry (this disallows 
shifts in the vev of S to absorb the linear term). However, for Nf < Nc this symmetry is 
anomalous because of the dynamically generated ADS contribution, and the R-symmetry 
is broken to a discrete subgroup. 

One may wonder this model breaks SUSY globally. If this were to happen (SUSY is 
unbroken), we could integrate out the gauge degrees of freedom and use the low energy 
effective description with mesons: 



A 7 

W = hS(M 12 - /2 2 ) + mi M n + m 2 M 22 + — — - (31) 

aet{M) 



The F-term equations are: 



F s = h(M 12 - /i 2 ) 

W — m A 7 M 22 
i'Mn - mi - 

Fm 12 = hS + (32) 

j? A 7 Mi2 
FM ^ det(M) 2 

*Mu -m 2 - det[M)2 



These equations can be solved asymptotically with 



C _ A 7 M 2 i 
° hdet(M) 2 

M_ m,2&? 
11 — € 2 A T 
,2 



M i2 = ^ (33) 



M 22 = ^ 



where A is some finite mass scale and e — > is some dimensionless parameter. Note 
that in this limit all vevs but M\ 2 go to infinity. 
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In the D-flat directions the Kahler potential is given by: 

K = VMtM (34) 

Where this is valid as long as the theory is weakly coupled (i.e) mesonic vevs are larger 
than A. 

The second derivative of the Kahler potential goes to as e — »• 0. One can, for 
simplicity take mi = mi = m, diagonalize M with unitary transformations and do the 
computation explicitly (using Mathematica for example) and see that in the runaway 
direction the V F = K~ b 1 W a W b goes to 2m 2 /i 2 . 

This can be understood since as e becomes very small, the scale at which the gauge 
symmetry is broken becomes large and the theory is weakly coupled. The reason why the 
potential doesn't slope to for very large vevs of the meson field is easy to understand 
since in this regime the theory should be weakly coupled and the microscopic description 
should be good. Furthermore, in this regime, the instanton contribution is rather small 
compared with the tree-level superpotential, and there is an approximate R-symmetry 
that guarantees that SUSY is broken. 

The reason why we can't just expand the Kahler potential around some scale to get a 
"canonical" Kahler potential for the meson fields (normalized by the square of the scale 
we are expanding around), is that there is no scale around which we can do this, i.e. there 
is no scale around which the minimization of Vf will justify the fact that we've neglected 
higher order terms in the expansion of the Kahler potential (Vp will be a function that 
slopes to when the fields go to infinity). 

Witten Index 

One may wonder if the statement that this model breaks SUSY is in disagreement with 
the Witten index argument, but note that taking mi and m 2 to infinity (the mass matrix 
has maximal rank) allows us to integrate out the quark fields (m susy <C mi, 7712). The 
scalar field is then light and survives to low energy, doing this integration one sees that 
at low energies the effective superpotential is W = —hSjJ 2 + NcA 3 SYM . In this case the 
Det(M) = mim 2 , where M is the quark mass matrix (as in [23]), and A$ym does not 
depend on the vev of S. There is no mass term for S in the superpotential and because of 
the linear term, the model breaks SUSY. So the Witten index must vanish. 

Numerical Calculations 

Using Mathematica one can compute the eigenvalues of the matrix M^M and do the 
square root. The sign is determined by the physical requirement that the potential should 
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be bounded from bellow. It's straightforward to compute the Kahler metric and invert 
it to compute the Vf- One can then give values to the parameters and Mathematica 
provides different numerical algorithms that help finding the minimum of the potential. 

As a test example, we use mi = 1.2,rrt2 = 0.45, h = 0.8, A = 0.05 in units of /i, and 
g = 0.75. Using the mesonic degrees of freedom, we find that there is a minimum for: 



S = -23.36 
M n = 2.43 

Mi 2 = 0.16 (35) 
M 21 = 101.10 
M 22 = 6.49 

with Vf = 0.62/i 4 . However, this minimizes Vf along the directions where Vd = 0, 
and may not be the right procedure. Instead of improving our understanding of the 
Kahler potential around the classical D-flat "pseudo-moduli space" , as in [TQl E] , we 
follow directly to the microscopic description. 

We now minimize the potential using the microscopic description, and allowing for 
non-vanishing D-terms. The Vf potential is: 

\r !,2|j ~ 2 1 2 I I . A 7 (X-X0i - <t>-XXi) ,2 . |r C ~ , 1 A? (X-X0i - X-iXi) 12 

V F = h \<p.x - [i \ + I m i0i = — = — —\ + \hSXi + m\<?i = — = — —\ 

{<P4x-x - x-H-x) 2 w4x-x - x-H-x) 2 

Hi) ,2,1 ~ A 7 (0.0Xi - 4>-x4>i) 

— Z — I + \ m 2Xi ~ Z = 

■x-x4$-x) 2 {$4x-x - x-H-x) 2 ' 



,|l i , A 7 (0.0Xi - X-Hi) ,2,1 ~ A 7 (0.0Xi - ,2 (oa\ 



Where a sum over i is assumed. We chose some gauge fixing conditions: <p2 — 0, 
03 = 0, X3 = and minimized the potential with respect to the other fields, where the 
lower case index is a color index. Numerically we found that 03 = 0,X3 = 0. 

The non-vanishing vevs are given (in units of fi) by: 



s = - 


-1.51 x lO" 3 


4 


»i = 0.50 


1 = 


= 1.4 x 10~ 3 


02 = 


1.36 x 10" 4 


Xi = 


1.66 x 10^ 3 


X2 = 


-2.28 x 10~ E 



(37) 



Xi = 0.35 
X2 = -0.67 

And Vf + Vd = 0.58/i 4 , which is lower than the previous result. This vacuum has 
Vp = 0.54/i 4 and Vd = 0.04/i 4 , and so both D and F-terms are non-zero. This means 
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that by moving away from the "D-flat" directions we were able to find a global minimum 
where SUSY is broken, and because the scale at which the gauge group is broken is much 
larger than A, the theory is weakly coupled. 

We note that the scale at wich gauge symmetry is broken is only one order of magnitude 
above the strong coupling regime, so quantum corrections are not negligible. One can 
change the parameters (for example, decreasing g will increase the vevs of the quarks), 
however, this does not change the fact that SUSY is broken. As a consequence of this 
fact, this model should be taken in the sense of an existence proof more than anything 
else. 

To summarize what we have done: we have took a model that had an R-symmetry 
and wrote the most general superpotential consistent with this. We noticed that for 
Nf < N c there are non-perturbative corrections to the potential that in practice make 
this symmetry anomalous. Despite having no continuous R-symmetry the model breaks 
SUSY because the superpotential is not generic and, for some range of the parameters, 
the global vacuum of the theory has a combined F and D-term SUSY breaking vacuum. 

The way this work makes contact with the work of [191 EQ] is that we can look at 
the gauge singlet as a Lagrange multiplier that deforms the "classical moduli space" of 
SQCD. On the other hand, the singlet also changes the behavior of the potential at infinity 
allowing the Witten index to change from the SQCD case. One can readily see that in 
this model it has to vanish (by looking at the case where the quark masses are much 
larger than the SUSY breaking scale). 

A detailed study of the phenomenology and inclusion of messengers in this model is 
postponed for later. Some recent work in semi-direct mediation was done in [H] using 
very similar models to the one we presented here, except that they use it for the messenger 
sector. In this work they were able to retrofit the model to generate gaugino masses at 
leading order in F/M. A different approach to mediating the SUSY breaking effects to 
the MSSM was done in |12j (the messengers coupling to the hidden sector is an irrelevant 
operator suppressed by a high mass). 



4 Getting D-terms at loop level, the runaway case 

4.1 The model 

We will follow refs.JH El [6]. In the magnetic description of the theory, the baryon 
deformation consists in having an ISS model with magnetic gauge group N = 2 and 
Nf = 7, where we've added to the superpotential a quadratic term in the quark field: 
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W = h&utpi.tpi - nl^a) + me i3 e a \y h (38) 

The last term is a baryon of SU(2) and is a relevant deformation of the low energy 
magnetic description. The quark fields have a charge that is not or 2, allowing for 
spontaneous R-symmetry breaking. 

The masses /x^ break flavour symmetry from SU(7) down to SU(2) x SU(5), and can 
be written as: 



, _ / p% \ 
^ ~ { fi% J 

with fi 2 > (i 2 to keep the SU(5) symmetry unbroken. We can write the fields as: 





SU(2) f 


517(5), 






M<# + 1 □ ^ 


(1 □ \ 
V □ Adj + 1 J 


2 




D (?) 


(a 


2 
3 




(?) 


(a 


2 
3 



Classically there's a runaway direction parameterized by: 



,2 . 



when £ goes to oo. 

Motivated by these constrains, we will minimize the potential along the directions: 

And check that the pseudo-moduli have positive masses in the minimum we find. 

The mechanism through which SUSY breaking is communicated to the SSM is direct 
gauge mediation: we gauge the unbroken SU (5) of flavour and the fields charged under 
this symmetry will act as messengers. These will be of two types: fundamentals(p,p,Z 
and Z) and the adjoint X. The fundamentals will generate masses for the gauginos at one 
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loop and masses for the scalars at two loops in SSM, while the adjoints will contribute 
to gaugino masses through an effective one loop diagram [6]. For this particular model 
the contribution of the adjoints to gaugino masses is negligible when compared to the 
contribution of the fundamentals (even when we take into account the effects of D-terms 
and their derivatives). (It was recently noted that in the case where the gauge group is 
SO(N) there is no runaway and this deformation probably does not play a role in getting 
bigger gaugino masses). For simplicity, in the results presented, we will only show the 
contribution of the fundamentals. 

We define reduced masses mi/2 and mo for the gauginos and the scalars: 

m XA (fi) := ^m 1/2 

Where C A and Sa are the Casimir and Dynkin index of the gauge group. And m^ A (/i) 
and m,f(fj,) can be computed from the usual diagrams. 



4.2 The computation of the effective potential 

The Coleman- Weinberg potential is given by the following expression: 



(Tr(mtJog(-fL) _ 2 Tr(m 4 f %(-^) + 3TrK%(-f)) (43) 
2 ^uv A c/v iV UV 



The mass matrices are given by: 



_ / W ab W bc + D aa D a c + D° a D a W abc W b + D aa D ac 
~ I W abc W b + D a a D a c W ab W bc + D a a D ac + D™D a ' ' 



where W c = |^ is the derivative of the superpotential with respect to the scalar com- 
ponent of the superfield $ c , D a = gz a T h aa z b and D% is the derivative of the D-term with 
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respect to $ a . And, as before, indices are raised and lowered by complex conjugation. At 
one loop the effective potential is V e ff = Vq + Vcw ■ 

The difference of the computations presented below is the inclusion of D-term and their 
derivatives in the calculation of the Coleman- Weinberg potential. In other words, we take 
into account not only the dependence of the effective potential on the superpotential 
parameters but also the ISS gauge coupling. 

As we've said before, the reason why this is important is that the masses of the 
particles in the Higgsed vector multiplet will not be supersymmetric (even in the absence 
of D-terms), since they couple to some F-terms that are not gauge singlets. This effect 
will then feed into the Coleman- Weinberg and the minimization of the potential. The 
parameter that controls the relative size of these contributions to the ones coming from 
chiral fields (that were not Higgsed) is k = (f) . This is because contributions to the 
Coleman- Weinberg scale with F x m 2 , where F are the F-terms that break SUSY, and m 
is the mass of the particle; for chiral fields this is h 4 and for fields in the vector multiplet 
this is h 2 g 2 , and the ratio is the parameter k defined before. If k is large (but both h 
and g are small), contributions coming from the Coleman- Weinberg could change the tree 
level constraints by a sufficient amount to allow the generation of larger gaugino masses. 

When computing the one loop correction one needs to know the masses of all particles. 
In this particular model, because we have a runaway direction, this presents a problem. 
We have some global symmetries that are spontaneously broken. For each of the broken 
symmetry generators there should be a goldstone boson. However since we are not at the 
tree level minimum of the potential, these modes will not be massless. This only happens 
after the stabilization of the potential by quantum corrections. These will shift the mass 
of the modes corresponding to broken global symmetry generators by a small amount. 
What makes this relevant for these modes is that since this shift can be positive, their 
tree level masses can be negative. This is just a reflection that only the full quantum 
corrected potential has a well defined vacuum. In view of this we will not try to ensure 
the positivity of the masses for these modes (which we will call Goldstone modes anyway), 
and estimate the error by looking at the minimum of the following functions: 

V cw = V cw (Massive modes) 

Vcw = V CW (M as sive modes) ± V cw {Gold stone modes) 

This is similar to including the one-loop corrected mass for Goldstone modes and tree-level 
masses for all other fields when computing the Coleman- Weinberg potential. 

By looking at several cases with different parameter choices, we've seen that usually 
the error introduced in this approximation is small, and probably of the same order 
of two loop corrections. We show two representative cases (the diagonalization of the 
mass matrices was done numerically and the error is omitted whenever it's smaller than 
numerical sensitivity): 
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k 


V 


K 


z 


X 





-0.1601 


0.5674 ± 1 x 10~ 7 


7.0403 ± 3 x 10~ 6 


-0.3157 ±4 x 10~ 7 


0.36 


-0.1602 


0.5692 


7.0570 


-0.3177 


1.0 


-0.1566 ±2 x 10~ 7 


0.5671 ± 1 x 10" 7 


7.1914 ± 1 x 10" 6 


-0.3155 ±5 x 10" 6 



Table 1: Values of the vev's for the baryon deformation using h = 1, m — — 2 and 
A = 1, A:= (|) 2 . 

4.3 Results 

We now study what happens when we switch on the gauge coupling: 

As one could expect, and can see in tables 1-3 the position at which the vacuum 
is stabilized depends on the gauge coupling. The violation of the classical constraints 
increase as we increase k (there is a decrease at the beginning because h and g contributions 
to the violation of the tree level constraints have opposite signs). There is, however a 
maximum value for k, after which the potential is not stabilized at one loop, for the point 
that we show this is k ~ 4 . The reason why this happens is that the F-terms that 
exist along the runaway cause a mixing in the fermionic mass matrix between the N 2 — 1 
gauginos of the hidden sector and the fermions. This is proportional to g 2 h 2 and comes 
with a factor of -2 in the Coleman- Weinberg. In the scalar mass matrix, there is a mixing 
of the N 2 — 1 "Higgsed" particles proportional to g 2 h 2 and this comes with a factor of 1. 
There is also the usual mixing due to F-terms and this is proportional to h A . So, there 
are directions in field space where the Coleman- Weinberg behaves as h A (l — ak). If a is 
large enough the potential will not be stabilized. 



k 


V 






X 





-0.1634 


0.3617 


11.0533 


-0.3258 


0.44 


-0.1635 


0.3622 


11.0660 


-0.3268 


1.78 


-0.1578 


0.3581 


11.3329 


-0.3194 


3.36 


-0.1259 


0.3229 


12.8147 


-0.2600 



Table 2: Values of the vev's for the baryon deformation using h = 0.6, m = 1.5, /i — 2 
and pL—1. 

We now turn to the computation of the ratio of gaugino vs scalar masses. Before 
going to the numerical calculations we can understand what will happen by noting that in 
leading order of F/M, gaugino masses vanish when the classical constraints are satisfied. 
So we expect that they increase when the violations of these constraints increase. By 
looking at table 3, we see that we expect a decrease of the ratio of gaugino vs scalar mass 
for small k and then an increase. This is indeed what happens: 
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k 


Vc-V 




mi 

■>■ 





3.2 x 10~ 5 


1.69 x 1(T 4 


3.85 x 10" 4 


0.44 


4.55 x 10" 4 


-7.7 x 10" 4 


4.02 x 10" 4 


1.78 


2.36 x 10~ 3 


-5.16 x 10- 3 


5.29 x 10" 4 


3.36 


4.23 x 10~ 3 


-1.08 x 10~ 2 


8.00 x 10" 4 



Table 3: Evolution of the violation of constraints with k and soft gaugino masses. 



k 


m 


mi /mn 

2 





0.402 


9.58 x 10" 4 


0.44 


0.402 


1.00 x 10~ 3 


1.78 


0.394 


1.34 x 10~ 3 


3.36 


0.356 


2.23 x 10~ 3 



Table 4: Change of the ratio of gaugino to scalar reduced masses with k, h = 0.6, m = 1.5, 
[i = 2 and p, — 1 

One could hope to reduce the hierarchy qualitatively between gaugino and scalar 
masses when k > 1, but even though there is a decrease it is not enough: unfortunately, 
in the promising regime of large k (but both h and g are small), we stop being able to 
stabilize the potential with one loop corrections. 

In these cases, even though there are effects coming from mixing derivatives of D-terms 
and F-terms, there are no non-vanishing D-terms because of an unbroken SU(2) V . Under 
this global symmetry the quarks transforms as fundamental x anti-fundamental, i.e. the 
only vev it can get is for the singlet component, and this must be proportional to the 
identity matrix. This symmetry could have been spontaneously broken but, numerically, 
we didn't find this to happen. 

We generate non- vanishing D-terms by breaking the SU (2) f global symmetry by mass 
terms, which also breaks the SU(2) V . These contribute to the tree level potential, so one 
might thing that they will not be generated. The reason why they do appear, is because 
there is a decrease in the Coleman- Weinberg that compensates for this increase. The mass 
matrix is now: 

//i? \ 
/4 = A (46) 
V jl% J 

The vev is parameterized as before except that now 0, and Y are 2x2 matrices 
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Figure 1: Impact of the gauge coupling on the ratio — where m A is the reduced gaugino 
mass and m s is the reduced scalar mass, for m = 0.3 in red, m — 1 in blue, m = 1.5 in 
green, m = 2 in purple and m = 2.5 in black. 



with different diagonal values: 



«i \ 7 / 6 

« 2 7 ' 9 V & 



(47) 



We have checked that non-diagonal vev's for 0, and Y are not generated after the 
inclusion of the Coleman-Weinberg potential, hence we used this simpler form for the 
minimization. 



With this ansatz for the vacuum only the third component of the D-terms D 3 can 
be non-zero. Also here we don't see a qualitative difference and gaugino masses remain 
small. 



k 




k 2 


6 


6 


Vi 


V2 


X 





0.66 


0.30 


13.63 


13.29 


-1.01 x lO" 1 


-2.27 x lO" 1 


-3.37 x 10- 1 


0.735 


0.67 


0.28 


14.60 


14.58 


-0.99 x 10- 1 


-2.24 x lO" 1 


-4.71 x 10~ 2 


1.306 


0.64 


0.28 


14.31 


14.29 


-0.91 x 10- 1 


-2.04 x 10- 1 


-4.08 x 10- 2 



Table 5: Change of the vacuum position with k, for /ii = 3.0, /i 2 = 2.0, jl — 1, m — 1.6, 
h = 0.7 
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k 


D 3 /F 


mi Irrin 

2 





0.0 


1.20 x 10~ 3 


0.735 


-3.2 x 10~ 3 


1.21 x 10" 3 


1.306 


-4.1 x 10~ 3 


1.61 x 10~ 3 



Table 6: Change in gaugino mass and D-terms with k. The D-term was rescaled by the 
SUSY breaking scale. 

5 Conclusions 

The motivation for this work was to try to evade the results of [7j, where they show 
that when the superpotential is renormalizable and SUSY is broken through F-terms, 
metastability is unavoidable if one wants to use gauge mediation to generate large gaugino 
masses in the MSSM. 

We investigated the possibility of having models that break SUSY through F and 
D-terms (very much in the style of pilj [TTj). It was seen that this does not happen if 
the superpotential is renormalizable and the Kahler potential is canonical. This also will 
not happen if the Kahler potential is canonical and the non-vanishing F-terms are gauge 
invariant. (Metastable vacua with non-vanishing F and D-terms are possible, but care 
needs to be taken regarding flat directions since often these will not be stabilized by one 
loop corrections.) 

We showed that when these condition are not met one can have models that break 
SUSY in the global vacuum of the effective description with both non-vanishing F and 
D-terms: SQCD+singlet and baryon deformed ISS. 

In the first case, and by changing the parameters, one goes from a pure F-term SUSY 
breaking to a combined F and D-term SUSY breaking model. This is not in contradiction 
with the Witten index argument since the presence of the singlet changes the potential at 
infinity. In particular, one can choose a region of parameter space where we can integrate 
out the quarks and at low energies one has a gauge singlet with a classically flat (but 
non-zero) potential. 

For our particular choice of parameters, gauge symmetry was broken only one order 
of magnitude above the dynamical scale of the gauge theory, so quantum corrections to 
the potential are not negligible. This could be improved by changing the parameters (eg. 
reducing the gauge coupling increases the size of some of the quark vevs). Due to the 
existence of several tree-level minima with similar energies, study of Coleman- Weinberg 
corrections is very likelly to be important. This is left for future work. 

This model does not include messenger fields however. If one allows for metastability, 
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one can couple a vector-like pair of hidden sector singlets and SU(5) fundamentals to 
the gauge singlet and have ordinary gauge mediation (see [26] for a review). It may be 
possible to do things in such a way that it makes the discrete R-symmetry[27j and D-terms 
present in the model useful. 

This model is another example that combined F and D-term SUSY breaking can be 
achieved (in a theory with vector-like matter), but more work is required to build a 
phenomenologically viable model. 

The second case is fairly similar to the first, but the vev's of D-terms are loop supressed 
with respect to the F-terms. One could hope that by increasing the ISS gauge coupling 
with respect to the ISS parameter h, the departure from the conditions in [7] would be 
sufficient to generate larger gaugino masses. Unfortunately, in this case, and as in one 
of the models presented in [21], there is a maximum value of g for which the runaway is 
stabilized. For this model, similar gaugino and scalar masses will not be achieved without 
some tuning. 
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